Long-Term Stability of an Area-Reversible Atom-Interferometer Sagnac Gyroscope 
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We report on a study of the long-term stability and absolute accuracy of an atom interferometer 
gyroscope. This study included the implementation of an electro-optical technique to reverse the 
vector area of the interferometer for reduced systematics and a careful study of systematic phase 
shifts. Our data strongly suggests that drifts less than 96 pdeg/hr are possible after empirically 
removing shifts due to measured changes in temperature, laser intensity, and several other experi- 
mental parameters. 
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A Sagnac gyroscope measures rotation through an in- 
duced phase shift, known as the Sagnac phase, between 
the two arms of an interferometer. In recent years Sagnac 
gyroscopes have demonstrated extraordinary sensitivity 
lily. L| . In addition to sensitivity, however, many appli- 
cations, including inertial navigation and tests of general 
relativity, require an instrument with high accuracy or 
excellent long-term stability 4]. In this paper we re- 
port on a study of the long-term performance of an atom 
interferometer gyroscope. As part of this study we im- 
plemented a method to periodically reverse the vector 
area of the interferometer loop, a technique which has 
already proven useful in atom interferometer accelerome- 
ters [3 • This reversal causes the sign of the Sagnac phase 
to change, allowing systematics which do not change sign 
with the reversal to be cancelled. 

The apparatus has been described previously Q, 0. 
Experimental details are summarized here. A pair of 
counter-propagating cesium atomic beams run the length 
of the apparatus. The atoms are generated thermally 
and have a characteristic velocity of 220 m/s (determined 
from the Sagnac phase induced by Earth rotation and by 
time-of- flight measurements). The atoms are collimated 
by narrow tubes, transversely laser cooled, and then op- 
tically pumped into the 6 2 S!/ 2 F = 3 hyperfine ground 
state. They then interact with three pairs of counter- 
propagating laser beams traveling perpendicular to the 
atomic beams. The three interaction regions are spaced 
by 0.968 m, for a total length from the center of the 
first to the center of the last beam pair of 1.936 m. The 
counter-propagating laser beams drive stimulated Raman 
transitions between the F = 3, mp = and the F = 4, 
mp = cesium ground-state hyperfine levels. Laser po- 
larization and a magnetic bias field suppress transitions 
between other magnetic sub-levels. For the studies re- 
ported here the two Raman lasers were tuned about 850 
MHz below the 6 2 S 1/2 F = 3 and F = 4 to 6 2 P 3/2 F = 3 
transitions. This combination of parameters results in a 



Sagnac phase of 9.1 radians for a rotation rate equal to 
Earth rotation (f2g = 15 deg/hr). 

Each Raman transition involves the absorption of a 
photon from one laser beam and the stimulated emission 
of a photon into the counter-propagating beam, giving 
the atom two photon recoils of transverse momentum. 
As a result, the Raman lasers act as atom "mirrors" and 
"beamsplitters" to spatially split and then coherently re- 
combine the beam of atoms. In the language of Ramsey 
interferometry the laser fields form a7r/2 — it — ir/2 pulse 
sequence. The first ir/2 pulse places the atoms into an 
equal superposition of F = 3 and F = 4 states and gives 
the F = 4 component two photon recoils of momentum. 
Next the tt pulse swaps the internal energy and external 
momentum states of the two halves of the atom wave. 
This causes their paths to cross at the final 7r/2 pulse 
which mixes them and causes them to interfere. The 
atoms are then probed by resonant fluorescence to deter- 
mine the final population in the F = 4 state. From this 
measurement the interferometer phase, and thereby the 
rotation rate of the apparatus, is determined. 

To implement the area reversal we modified our previ- 
ous method of generating the two Raman laser frequen- 
cies 0- In our current scheme, the +1 and -1 diffraction 
orders of a master laser passing through a high frequency 
acousto-optic modulator (AOM) are used to inject two 
slave diode lasers. The AOM's frequency is set such that 
the slave lasers differ in frequency by the cesium hyper- 
fine ground state splitting plus or minus a small (5 MHz) 
detuning. The 5 MHz detuning from the two-photon 
resonance prevents the co-propagating light frequencies 
from driving recoil- free Raman transitions. The beams 
from the two slave lasers are coupled into a common fiber 
which delivers the light to the interferometer. The beam 
is then split into three "Raman" beams: two of equal 
intensity plus a third beam at twice this intensity. These 
beams are used to generate the tt and 7r/2 pulses. 

The sign of the 5 MHz detuning determines the di- 
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rection of the interferometer's vector area. Each Ra- 
man beam passes through the atomic beam in one of 
the three interaction regions. After passing once through 
the atomic beam, each Raman beam is sent through a 
pair of AOMs and then retro-reflected back through the 
AOMs and the atomic beam. The double pass through 
the AOMs results in a frequency upshift of 5 MHz. This 
causes one of the incoming light frequencies to form a 
two-photon resonance with one of the retro-reflected fre- 
quencies. Changing the sign of the 5 MHz detuning 
switches which frequency components of the direct and 
retro-reflected light are involved in the transition. This 
causes the atoms to recoil in the opposite direction, re- 
versing the interferometer's area. This electro-optic area 
reversal method is extremely clean, producing a pre- 
cise 180° reversal and exceptional cancellation of numer- 
ous systematic errors which change sign relative to the 
Sagnac phase with the reversal. 

All of the data we present was taken utilizing the 
counter-propagating atomic beam and electro-optically 
applied rotation bias methods described in pj, allow- 
ing us to cancel linear acceleration shifts and the Sagnac 
phase due to the Earth's rotation. The ability to bias our 
instrument near zero rotation is an uncommon feature 
in Sagnac interferometers. Laser gyros typically require 
a rotation bias to prevent backscatter-related frequency 
pulling and mode locking ||. Working near zero bias 
makes the interferometer insensitive to many sources of 
instability and increases the fringe contrast for improved 
signal-to-noise ratios. In this mode of operation the in- 
strument's rotation rate can be obtained by subtracting 
the electro-optically applied rotation bias from the inter- 
ferometrically measured rotation rate. 

To study long-term stability we took several sets of 
data, each of which span a period of about two and a half 
days. A typical data set is shown in Fig.Q] Although this 
data set exhibits more long-term drift than some of our 
other sets, it demonstrates the cancellation of system- 
atic drifts that can result from area reversal. This figure 
presents raw data, corrected only by the removal of a cali- 
brated, independently-measured phase-shift due to drifts 
in Raman laser intensity. The strong anti-correlation of 
the two configurations seen in this figure, and the much 
lower drift in the average phase relative to either con- 
figuration, demonstrates the utility of area reversal for 
increased stability. 

Figure ^ shows high-frequency noise on top of slow 
drifts. The high-frequency noise is due to mechanical 
noise in the building. For measurements of long-term 
stability we placed the interferometer directly on the lab 
floor to reduce measurement errors due to slow drift in 
the position of the device. Although this results in ef- 
ficient coupling of our apparatus to the high-frequency 
mechanical noise present in the room, for these studies 
we were not concerned with high-frequency noise which 
quickly averages away. This noise has been well char- 
acterized with high-speed measurements of the Sagnac 
phase of our instrument. The rotational noise, shown in 
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FIG. 1: Variation of measured phase over a 60 hour period. 
The upper and lower traces are the interferometer phase in 
the area reversed and non-reversed configurations (± a con- 
stant such that they don't overlap) . The interferometer phase 
of the reversed-area trace has been multiplied by -1 such that 
a positive phase change indicates a positive rotation change in 
both traces. The middle trace is the average of the other two 
traces. Rapid drifts at the beginning of data sets, evident in 
the upper and lower traces, were very common, possibly due 
to the change in temperature after the experimenters left the 
room. Area reversal does an extremely good job of remov- 
ing this drift. For comparison, the Sagnac phase for Earth 
rotation rate is 9.1 radians. 
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FIG. 2: The power spectral density of the random angular 
motion. This plot shows the typical rotational noise present 
in our laboratory. 



Fig. El hits a baseline angle random walk (ARW) of a few 
times 10 -6 deg/hr 1 ' 2 at frequencies from 2 to 7 Hz. 

Several of our data sets showed considerably less low- 
frequency drift than the data shown in Fig^due to lower 
environmental noise in the lab during these runs. One of 
these data sets is shown in Fig.^a). The Allan variance 
of the raw data in this figure, shown in Fig. Etc), bottoms 
out at a measurement time of about 18 minutes. Using 
the methods described in we calculate a bias stability 
of 560 ^deg/hr for this data set. This is considerably 
smaller than what is found from the data in Fig. ^ which 
bottoms out at 1.1 mdeg/hr at a measurement time of 
about an hour. 

The stability of most gyroscopes is improved consid- 
erably by removing known correlations in the rotation 
signal with other measured parameters such as the tem- 
perature of the apparatus. In each experiment we mea- 
sured ~ 30 auxiliary parameters, including laser pointing 
and intensity and temperatures around the apparatus. 
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FIG. 3: The effect of removing correlations. The black line in 
(a) is the Sagnac phase from a particularly quiet set of data 
plotted vs. time. This data represents the average rotation 
phase for the area reversed and non-reversed configuration. 
The thick grey line shows the correlation of the Sagnac phase 
drift with 9 other measured parameters. The line represents 
the sum of the parameters each multiplied by a fixed constant. 
The plot in (b) is the same Sagnac phase data after the data is 
corrected by subtracting these linear correlations. Plot (c) is 
the Allan deviation of the Sagnac phase. The dashed line rep- 
resents the deviation of the uncorrected data. The solid line 
represents the corrected data. The dotted line results from 
correcting the Sagnac phase by removing linear correlations 
with just three parameters. 



In quiet data sets the long-term drift shows strong corre- 
lation with several of these parameters. A simple linear 
correction, in which a subset of these quantities were each 
multiplied by a constant and subtracted from the Sagnac 
phase at every point in time, considerably reduced drift. 
The effect is shown in Fig. Ota) and (b) . The Allan de- 
viation is plotted in Fig.^c). 

The Allan deviation in Fig. |3{c) shows that the cor- 
rected signal's long-term stability is limited not by sys- 
tematic error but by random noise, even for the maxi- 
mum measurement time of 4.7 hours. Because the last 
few points in the Allan variance involve only a statisti- 
cally small number of clusters of data [9|, rather than 
use the minimum Allan deviation of 22 /xrad, we calcu- 
lated a more conservative estimate of the bias stability of 
the corrected data by extrapolating out to 4.7 hours the 
i -1 / 2 behavior present in Fig. IHIc) up to measurement 
times of about 1 hour. This yields a minimum devia- 
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FIG. 4: Cancellation of systematics by area reversal. The 
circles/squares represent the reversed/non-reversed configu- 
ration. The measured interferometer phase is plotted verses 
(a) the applied magnetic bias field, and (b) the center pulse 
displacement. The phase shift with magnetic field in (a) is 
parabolic, as expected. The offset of the apex of each parabola 
from zero bias indicates the presence of stray magnetic fields. 
The filled and unfilled marks in (b) represent measurements 
on atomic beams with two different average horizontal trans- 
verse velocities. Averaging the reversed and non-reversed 
phase cancels the magnetic bias systematic and most of the 
center beam displacement systematic. 

tion of 88 ^trad corresponding to a bias stability of 96 
/j,deg/hr. The character of the Allan variance strongly 
suggests that the inherent long-term drift of the appara- 
tus could be less than this and implies a scale factor drift 
of less than 7 parts per million. 

In addition to stability tests, we have also re-examined 
systematics in the context of area reversal 0, |(| . Below 
we discuss only the dominant sources. A detailed discus- 
sion of the theory can be found in ^(J • 

One significant error results from a ~ 0.5% difference 
in the magnetic bias field between the two halves of the 
apparatus. The resulting quadratic Zeeman shift system- 
atic is shown in Fig. 0[a). Area reversal nearly cancels 
this shift, leaving a residual shift of less than 0.1% of CIe- 

Another significant systematic phase shift occurs if the 
7T pulse is not exactly centered between the 7r/2 pulses. 
Assuming that gravity is nearly perpendicular to the 
plane defined by the interferometer and that the atomic 
beam is nearly perpendicular to the Raman beams, to 
lowest order this interferometer phase shift is given by 

= / hk eff 2 + 2k eff ■ v + ALk eff ■ (fl x v) | 
\ mv v v 2 

2Lk eff -g \ A 
v 2 J 

Here k e f / is the effective two-photon wave vector for the 
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Raman beams (denned such that the two-photon momen- 
tum kick the atom receives is given by hk e ff) and L is 
the distance between pulses. The parameters v and m are 
the effective velocity of the atomic beam and the atomic 
mass, and fi and g are the rotation rate and acceleration 
(gravity) present on the apparatus. The displacement of 
the 7r pulse from the center of the interferometer is rep- 
resented by A. The first term in Eq. ^ is cancelled by 
area reversal. Given our ability to center the it pulse, this 
term can account for most of the difference between the 
reversed and non-reversed phase near zero field in Fig. 
0Ja). The second term (a combination of A and a non- 
zero transverse atom velocity), and the last two terms 
(resulting from a combination of incrtial forces and un- 
balanced lengths), are not cancelled by area reversal. We 
studied these systematics by varying the pulse spacing. 
From the results, shown in Fig. EJb), we conclude that 
they contribute an error of less than 1% of 

Systematics due to horizontal Raman beam misalign- 
ment should be negligible; we use observed Doppler shifts 
in individual interaction regions for precision horizontal 
alignment. Vertical misalignments, however, may con- 
tribute significant shifts. Tilting the plane of the interfer- 
ometer introduces a phase shift due to gravity. This shift 
is cancelled using counter-propagating atom beams 0,03]. 
Deviation of the ir/2 beams from the plane defined by the 
atomic beam and the it beam results in a gravitational 
shift which is not cancelled by counter-propagating atom 
beams or with area reversal. A related shift of the same 
magnitude results from a combination of vertical Raman 
beam misalignment and a non-zero vertical component of 
an atom's velocity. Intentionally walking the vertical Ra- 
man beam alignment we could only change the measured 
rotation rate by 7% of £Ie before signal contrast was re- 
duced significantly. When aligned for maximal contrast 
this systematic should be below 2% of flE- 

Drift in Raman laser intensities introduce additional 
phase shifts. The two Raman laser intensities are cho- 
sen to minimize the differential ac-Stark shift between 
the F = 3 and F — 4 hyperfine ground-states. Be- 
cause both interferometer paths experience nearly equal 



ac-Stark shifts, and since Stark shifts can be cancelled 
using counter-propagating atomic beams, Stark effects 
should not be significant 0, 0] . But because the pulse 
area received by an atom with a given velocity depends 
on the intensities of the Raman beams, the laser intensity 
determines which velocities result in highest fringe con- 
trast and contribute most to the interferometer signal. 
Since the Sagnac phase at a given rotation rate depends 
on velocity, we still observe changes in phase when the 
intensities of the Raman beams are varied. 

Area reversal reduces intensity-dependent shifts by an 
order of magnitude. Incomplete cancellation is due in 
part to an imperfect implementation of area reversal; a 
difference in the detuning from the single photon res- 
onance for the two configurations and imperfect AOM 
efficiency in the retro-reflections results in slightly differ- 
ent pulse areas and ac-Stark shifts for the non-reversed 
and reversed configurations. Although it is possible to 
improve the implementation, these shifts are not a domi- 
nant source of instability. The difference in Sagnac phase 
for the two areas at a given rotation rate is well below 1%. 
And since we typically add a rotation bias to keep the 
Sagnac phase small, effects related to the atomic beam 
velocity which scale the Sagnac phase without adding an 
independent bias phase are greatly reduced. 

In conclusion, we have implemented a method of re- 
versing the vector area of the interferometer to reduce 
systematic drifts and have completed a study of the long- 
term behavior of our gyroscope. Removing correlations 
with auxiliary parameters reveals that the intrinsic sta- 
bility of our gyroscope is probably less than 96 ^deg/hr 
(this upper value limited by statistical white noise) with 
a scale factor stability better than 7 parts per million. 
We also completed an analysis of the dominant system- 
atics present in our apparatus and determined that the 
absolute accuracy is better than a few percent of Earth 
rotation. 
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